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^ Abstract 
i 

l A mathematical model of plastic deformation in face-centered cubic (FCC) materials based on a balance model taking into account 

> fundamental properties of deformation defects of a crystal lattice was developed. This model is based on a system of ordinary 

differential equations (ODE) accounting for various mechanisms of generation and annihilation of deformation defects for different 
external conditions. In-house developed software, SPFCC (Slip Plasticity of Face-Centered Cubic), was employed to solve the 
system of ordinary differential equations. The implemented code solves efficiently the stiff ODE system and provides a user- 
friendly interface for investigation of various features of plastic deformation in FCC materials. 

Simulation of plastic deformation in the FCC metals was performed for the case of constant strain rate. The modelling results 
were validated by comparing experimental data and simulation results (stress-strain curves) and good agreement was obtained. 
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1. Introduction 

<N 

The role and possibilities of mathematical modelling and computer simulation in the analysis of plastic defor- 
mation are of great importance. The mathematical modelling allows to reveal the role of various factors influencing 
processes of plastic deformation such as mechanisms responsible for strain hardening and evolution of deformation 
defects of crystalline materials [£T3-[3j . Furthermore, the mathematical modelling is increasingly utilised for planning 
of experimental studies. 

The mathematical modelling of plastic deformation in crystals is based on mathematical relationship taking into 
account fundamental properties of deformation defects of crystal lattice and represents a fundamental part of crystal 
• i-H plasticity models ll4ll6 HT0ll . The models usually consist of kinetic equation, which establishes a relationship between 

the current microstructure and the material's mechanical response. This relationship^ written in general form as 
ordinary differential equations (ODE). Various models reported in literature ll5l l6l l8l flTT - TP?) differ from each other in 
way the set of deformation defects and the mechanisms of their generation and annihilation are considered. 

Typically mathematical models of plastic deformation include the equations of deformation defects accounting 
for elementary processes and mechanisms of plasticity 131161171 1131 . Specifically, the choice of deformation defects, 
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defining the regularities of plastic deformation and the mechanisms of their generation and annihilation and mutual 
transformation, determines the adequacy and possibilities of a mathematical model. 

Crystallographic shear zone is selected as a basic structural element with which regard mechanisms of slip plas- 
ticity are considered in mathematical model of plastic deformation by slip @ [TO] [14l Q3] . The description of the 
mechanisms of formation of shear zones formation is based on the fundamental physical and topological properties of 
defects responsible for plastic mass transfer. The fact that all parameters of the model either have physical or geomet- 
ric interpretation or can be defined from physical considerations, or their values can be limited to a certain interval of 
may be specified, is a characteristic feature of the model and represent a significant advantage over phenomenological 
models. 

Plastic deformation of slip primarily produces dislocations of various types and point defects. The processes of 
initiations of defects during formation of crystallographic shear zones are the mechanisms defining their generation. 

Mathematical model of plastic deformation in FCC metals |[T6l [171 consists of the balance equations of shear- 
producing dislocations (their density is denoted as p m ), dislocations in vacancy (p l d ) and interstitial (p' d ) dipole config- 
urations, interstitial atoms (concentration c,), monovacancies (c\ v ) and divacancies (c2,,), the equation defining strain 
rate, and the equation describing external effects to the deformed material. The specific ingredients of the model of 
deformation defects generation and annihilation are derived from the assumption that: 

1 . The defects generation processes occur during the formation of crystallographic shear zones and are related to 
the dynamic nature of their formation. 

2. Annihilation processes have predominantly diffusion nature and occur in defect medium created by a set of 
defects generated by the slip in a large number of shear zones. 

3. Deformation defect continuum is homogeneous continuum and contains the same number of defects, that the 
entire crystallographic shear zones together. 

The mechanisms of deformation defects generation during the formation of crystallographic shear zone in monocrys- 
tals of FCC metals are reviewed in section 2, and the mechanisms of defects annihilation are presented in section 3. 
Section 4 summarises the system of differential equations of the model of plastic deformation in FCC metals. Sec- 
tion 5 is devoted to the description of SPFCC software, aimed to automatize the investigation of regularities of plastic 
deformation by slip with the use of a mathematical model. Section 6 gives an illustrative example showing possibilities 
of SPFCC software. 

2. The intensity of deformation defects generation in the formation of crystallographic shear zone 

The intensity of shear-producing dislocations generation. The intensity of the generation of shear-producing 
dislocations during the formation of crystallographic shear zone was obtained in [18| based on the ideas of Inden- 
bom V.L., Orlov A.K., and Kuhlmann-Wilsdorf D. |[T] |T5] |T9l . If we assume that the process of deformation occurs 
in identical shear zones, the intensity of the generation of shear-producing dislocations with deformation can be de- 
scribed as follows |T8l : 

G m (p) = F/(bD r ), (1) 

where D r = (B, t) / (Gbp m ) is an average diameter of the shear zone, b is the magnitude of the Burgers vector, parameter 
F is determined by the geometry of the dislocation loops and their distribution in the shear zone |18|, B, is the 
parameter determined by the probability of formation of extended dislocation barriers, r is the applied stress, G is the 
shear modulus. The description of the parameters and variables of the mathematical model is summirised in TabJT| 
(see Appendix A). 

Dislocations in dipole configurations are generated during the formation of crystallographic shear zone. The main 
cause for the occurrence of these dislocations is the presence of immobile dislocation jogs. When moving, the jogs 
leave a chain of point defects in case the height of the jogs is equal to one interatomic distance. In cases where the 
height of the jogs is more than one interatomic distance a dipole is formed. The dislocation loops emitted by the source 
broaden rapidly, due to decrease in nonlocal curvature. This leads to the generation of point defects with increasing 
intensity, and orientation close to screw, and their entrapment after short period of time. 



Gdiip) = Gdvip) = (6ydld(p)b) 
2 



(2) 
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where yj is the size parameter of the particles, Uip) is the length of the screw dislocation run. 

The intensity of generation of point defects. Non-conservative movement of dislocations with jog leads to 
the generation of such point deformation defects as interstitial atoms, vacancies and divacancies. The intensity of 
generation of point defects can be calculated as 

G k {p) = 2qT dy „/G, (3) 

where k = {i, lv, 2v}, Td yn is excessive stress (the stress excess over the static resistance to dislocation motion), and 
parameter q characterizing local concentration of stress during shear zone formation. 

If the distribution of forest dislocations is assumed random, and the intensity of vacancies and interstitial atoms 
generation is considered identical, the relation for the rate of the vacancies and the interstitial atoms generation is as 
follows: 

G,(p) = GivGo) = qTdyJG. (4) 

Taking into account the fact that the activation energy of divacancies migration is lower than the activation energy 
of monovacancies migration, and that divacancies migrate from vacancies chain behind threshold on dislocations to 
crystal volume presumably faster, the following relations for the intensity of mono- and divacancies generation during 
the process of deformation were derived: 

Giv(p) = qT dyn /(6G), (5) 
G 2v (p) = 5qT dyn /(6G). (6) 



3. Mathematical modelling of deformation defects annihilation processes 

Dislocations annihilation. A number of different mechanisms of dislocations annihilation have been presented 
in the literature ||4] [6] [7] [10] [16] [T8 1 . The main mechanisms of dislocations annihilation are: a) climb of non-screw 
dislocations as a result of precipitation of point defects on extra planes of non-screw dislocations, and b) annihilation 
of screw dislocations during their cross-slip. 

The following mechanisms are considered for annihilation of dipole dislocation configurations: 

1. Decrease in the shoulder of vacancy dislocation dipoles until their annihilation with precipitation of interstitial 
atoms onto them. 

2. Decrease in the shoulder of interstitial dislocation dipoles until their annihilation with precipitation of monovacan- 
cies and divacancies onto them. 

3. Increase of the shoulder of vacancy dislocation dipoles until the loss of their stability with precipitation of mono- 
vacancies and divacancies onto them. 

4. Increase of the shoulder of interstitial dislocation dipoles until the loss of their stability with precipitation of 
interstitial atoms onto them. 

Shear-producing dislocations annihilation by climb. Annihilation processes of diffusive nature in a dislocation 
subsystem include both climb of non-screw dislocations, and their counter-movement in the slip planes. Dislocations 
annihilation by crystallographic slip can occur only in pairs of dislocation segments, which are located closely enough 
to one another and have opposite Burgers vectors. The distance between such dislocation segments is defined by 
critical radius of capture: 

Gb 2-v 

" 47TT/ 1 — v' 

where Tt is friction stress, v is Poisson's ratio. 

At high density of dislocations when the distance between dislocations is less than r a , it can be assumed that the 

— 1/2 

distance between dislocation segments able to annihilate is equal to p m 1 . Assuming that dipoles are distributed uni- 
formly along the length of the shoulder and the distance between the dipoles in interval [0,p,„^ 2 ], then the expression 
for the rate of the shear-forming dislocations annihilation as a result of precipitation of interstitial atoms onto them 
can be written as: 

A m (p,c k ) = 2b' 1 p m min(r a , p m l/2 )aj' m /Ai(p,c k ), (7) 
3 
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where A,(p, q) = QkCk/a is the rate of point defects annihilation, k - {i, 1 v, 2v), Qk is the diffusion coefficient, a is the 
strain rate, a>' m is the fraction of interstitial atoms which left the glide plane on which they are stored. 

Screw dislocations annihilation by cross-slip. We assume that the cross-slip of dislocations occurs only at 
temperatures above a certain temperature T cs and is athermal. In this case, almost all screw dislocations annihilate if 
they are stable under current applied stress r and are trapped in dipole configurations in which the distance between 
dislocations of the opposite sign does not exceed the distance r a = Gb/(2jTTf). According to this approach the 
annihilation of screw dislocations by cross-slip can be represented by the relation: 

A£(pm) = FGp m co s P as /(BT), (8) 

where P as is the probability of annihilation of screw dislocations, a> s is a fraction of screw dislocations. 

Dislocations annihilation in dipole configurations. The main mechanism of dislocation annihilation in dipole 
configuration is the trapping of interstitial atoms on the vacancy dipoles, which leads to a decrease of the shoulder of 
the vacancy dipoles, and precipitation of interstitial atoms on interstitial dipoles. 

Dislocations annihilation in vacancy dipole configurations. Suppose that the maximum shoulder of the dipole, 
after which it loses its stability, is equal to the critical radius of the capture r a . In this case the rate of the dislocations 
annihilation in vacancy dipole configurations can be determined as follows: 

A dv (p,cd = 2cj?A i (p,c i )/(br a ). (9) 

Here a>. is the fraction of interstitial atoms, which left to vacancy dipoles. 

The density of dislocations forming vacancy dipoles during precipitation of vacancies onto them does not change 
significantly until the length of dipole shoulder becomes higher than r a . In this instance, dipoles loose stability and 
their further behaviour become analogous to that of the shear-forming dislocations. The rate of change of the density of 
dislocations in vacancy dipole configurations with deformation caused by precipitation of vacancies can be described 
as follows: 

R dv (p,c lv ) = 2co d l lA lv (p,c lv )/(br a l (10) 

where u dv v is the fraction of vacancies, which left to vacancy dipoles. 

The rate of change of the dislocation density in vacancy dipoles configurations connected with precipitation of 
monovacancies onto them can be described analogically: 

R dv (p, c 2v ) = 2< l A 2l ,(p, c 2v )/(br a ), (11) 

where coi^ is the fraction of divacancies, which left to vacancy dipoles. 

Dislocations annihilation in interstitial dipole configurations. Taking into account the fact that precipitation 
of vacancies onto interstitial dislocation dipoles leads to the decrease of the shoulder of the interstitial dislocation 
dipoles, whereas precipitation of interstitial atoms leads to the increase in their shoulder, similar to the intensity of the 
annihilation of vacancy dipoles considered above Eq. d9j, we can write down the following relations for the intensity 
of the annihilation of interstitial dipole configurations in case of precipitation of mono- and divacancies onto them 
respectively: 

A di {p,c u ) = 2a>f v A lv (p,c lv )l{br a ), (12) 

A di (p, c 2v ) = 2<A 2l ,(p, c 2v )/(br a ), (13) 

where o)f', oj d ' v are the fractions of mono- and divacancies respectively, which left to interstitial dipoles. 

For the intensity of change of the dislocations density in interstitial dipole configurations, which is connected with 
precipitation of interstitial atoms onto them, the following relation is used: 

Rdp, cd = 2ojfAi(p, cdl(br a ), (14) 

where co dl is a fraction of interstitial atoms, which left to interstitial dipoles. 

Annihilation of deformation point defects. During the process of plastic deformation, concentration of point 
defects is reduced as a result of their spending for the sinks. Dislocations, crystal grain boundaries, and crystal 
surface can act as the sinks for point defects. In this work, the following sinks are considered in particular models 

4 
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of point deformation defects annihilation: 1) for interstitial atoms - non-screw dislocations, mono-, divacancies; 
2) for monovacancies - non-screw dislocations, interstitial atoms, monovacancies; 3) for divacancies - non-screw 
dislocations, interstitial atoms. The rates of change of the deformation defects concentration can be expressed as 
follows: 

A k {p,c k ) = Q k c k /a. (15) 

It should be noted that, when two vacancies in their random movement meet, divacancies are formed, and when a 
divacancy meets an interstitial atom, or an interstitial atom meets a divacancy, a monovacancy is formed. Therefore, 
relations of deformation point defects accumulation take into account an additive to monovacancies and divacancies 
accumulation: 

Riv(c 2v , ci) = (Qi + Q2v)CiC2v/a, (16) 

*2»(civ) = Qwc\ v la. (17) 

The cases of point defects formed as a result of the meeting of two divacancies or a mono- and a divacancy were 
not considered in this work. 



4. The system of differential equations of deformation defects balance 

The presented models of deformation defects generation and annihilation described in sections 2-3 provide the 
basis for the creation of model of plastic deformation by slip. This model can be used to analyse the influence of 
various mechanisms and processes for FCC metals under various conditions of deformation. 

The system of differential equations of the balance of deformation defects for FCC metals can be summarized as 
follows 0: 



dC\ v _ qTdyn _ {Q\yC\ v - {Qi + ggyjcggy) 

da 6G a 
dc 2v _ 5qT dyn _ (QivCiv - Qivc\ v ) 
da 6G a 

dCj _ qTdyn _ QjCj 

da G a 

dp^ _ 1 _ 2(<jf "A,-(p, cd + <:A h ,(p, ci„) + <A 2l ,(p, c 2 „)) 

da 6y d ld(p)b br a 
dp d _ 1 2(wfj A u .(p, civ) + 0)f v A 2v (p, c 2v ) + <4jA;(p, c,)) 

da 6y d ld(p)b br a 
dp m F FGp m u s P as 2p m min(r a ,p m 1/2 ) a> k m 



(18) 



z 



da bDr Bt b k =\[jtM Ak< ^ ,Ck) 

2(<A h ,(p, Ch ,) + gj d 2 v A 2v (p,c 2v ) +pfA i (p,c i )) 
br a 

The last equation of the system ( fl"8] l takes into account the relaxation contribution due to the rate of change of the 
dislocations density in dipole configurations for the intensity of shear-forming dislocations accumulation. 

The equations establishing the relationship between strain rate the applied stress as well as the defect condition of 
the material are taken from ||H[T8). I n this work the following equation for the strain rate of plastic deformation is 
selected EJ: 

0,2Gb 3 -(T-T a )A(T,p)b 2 " 



fl = aoexp^ j, (19) 

where r„ is an athermal component of stress, Mj,p) is the function of stress and density of dislocations, Oq is pre- 
exponential factor. 

The system of equations ( fT8] > should be added by an equation (or equations), describing the influence causing 
deformation. 
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5. SPFCC software description 

A computational tool, denoted SPFCC (Slip Plasticity of Face-Centered Cubic), has been developed. This soft- 
ware can assist researchers in application of the mathematical models for the computational investigation of plastic 
deformation in FCC materials. Although, a number of commercial tool addressing this problem are readily available, a 
user needs to be familiar with methods of solving ODE and requires programming language skills in order to consider 
a specific problem. 

SPFCC software makes it possible to automatically form and calculate a full or reduced suite of defects of a 
mathematical model with little user intervention to set up the input values, run the calculations, plot stress strain 
curves, and export the results [ 18"]. SPFCC enables direct simulations of strain hardening and evolution of deformation 
defects of crystalline materials. 

Presently, the following types of effects can be considered: deformation with a constant strain rate a = const, at 
constant stress r = const and constant load P = const in conditions of tension t = Toexp(a/k s ) and compression 
t = TQexp(-a/k s ), where k s is Sachs factor. 

Fig.[T]shows main window of the software SPFCC. Simulation results (stress-strain curves) for different external 
conditions are obtained in less than one minute and posses a relative accuracy of the solution 10~ 3 . Fig.[T]illustrates the 
time of simulation results for constant strain rate at different temperatures. It should be noted that the computational 
time for the case of constant stress (creep) increases and the results are obtained in ten minutes. All the simulations 
were performed on the Supermicro computing cluster based on Intel Xeon Processor 5570 (code named "Nehalem", 
http://www.cluster.tpu.ru). 

SPFCC software helps to analyse the processes of plastic deformation in FCC metals and dispersion-strengthened 
materials with a metal FCC matrix and noncoherent, undeformable particles of the second phase 0] . 

The user can either apply the full mathematical model, or create reduced models. The reduced model can take 
into account the deformation defects and mechanisms of their generation and annihilation specified by the user. At the 
beginning of a calculation, the starting suite of defects and initial value of parameters are selected from the database. 

It should be noted that systems of ODE balance of deformation defects for various deformating effects and mate- 
rials may differ in the number of equations and the mechanisms of defects generation and annihilation because user 
can to choose or disregard some mechanisms in the reduced models. In addition, the equations can be piecewise 
cross-linked type. This places some restriction on the choice and realization of a numerical method. 



The system of differential equations of the model ( 18 H 19 1 is numerically stiff lfT7ll20] . This is mainly due to the 



fact that processes of generation and annihilation of deformation defects are significantly differ in terms of intensity. 
In addition value of variables are changing in the interval of integration by orders of magnitude. 

In order to solve a stiff ODE system, the software package SPFCC employs the explicit linear multistep Adams 
method to obtain initial values and variable order multistep Gear method (in the form of backward differentiation 
formulas). The Gear method is stable at any integration stepsize, and therefore can choose an integration step guided 
only by considerations of accuracy, but not stability |20|. The original Gear method has been modified 1171 . due 
to some physical restrictions imposed on the variables of the equations system (for example, variables can not be 
negative). The analysis of the effectiveness of numerical algorithms for the case of plastic deformation of dispersion- 
strengthened materials with undeformable particles was given in publication lTP7l . 

The graphical user interface of the SPFCC software allows to select variables and mechanisms of the genera- 
tion and annihilation of dislocations (Fig. [TJ. The user can chose values of the mathematical model parameters, 
initial conditions, and calculations accuracy in a dialogue mode. Furthermore, two material FCC metal or dispersion- 
strengthened material can be selected by the user.I n this case, the user will automatically be offered to choose the 
values of the parameters of material from the database software. It should be noted that the database is compiled using 
the results of theoretical and experimental studies available in literature. The user can accept the default values or set 
its own parameters values. If a user enters an incorrect values the SPFCC software will generate an error message and 
will prompt user to change value or to set the default values. 

The user can perform a series of calculations with the help of SPFCC software. To run the calculation, the 
parameter of the mathematical model (stress or load, temperature, initial density of dislocations, etc.), needs to be 
chosen together with the interval and the step of this parameters. 

In the process of the calculations it is possible to visualize intermediate results of modelling in real time, which 
allows to control the progress of calculations. The modelling results can be presented graphically. The user must 

6 
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Figure 1. Main window of SPFCC software. 



tick the checkbox in the table "Modeling results" and click next on "Reports" button in the toolbar. Fig. [2] shows an 
example of the flow curves for various FCC metals (aluminum, copper, nickel) at the temperature T = 293K. 

Modelling results and the values of the model parameters used are saved in a local database automatically. In- 
formation regarding the progress of numerical solution of the task is also recorded in the database. This information 
allows to generate recommendations for further computer experiments. 

The SPFCC software provides possibility to compare modelling results for various applied effects in tables and 
graphical form. According to modelling results, dependencies between any variables as well as temperature and strain 
rate dependency can be represented graphically. 



6. Illustrative example 

As an example of the model described above we show modelling results for the case of constant strain rate and its 
comparison with experimentally measured data. In case of the deformation of crystal with a constant strain rate, the 
equation of applied influence has the following form: a = const. This is a transcendental equation, which allows to 
find the stress at any given time of the deformation. 



Using the above-mentioned equations of deformation defects balance (18i and the equation defining strain rate 
( 19 1, calculations were performed for FCC metals with the values typical for copper: b = 2.5 ■ 10 10 m, vq = 10 13 i _1 , 



a = 0.5, Tf = IMPa, v - 1/3, k = 1.38 • 10 23 J/K, a> s = 0.3. The initial density of shear-forming dislocations were 
chosen to be 10 12 m~ 2 , and the initial density of dislocations in vacancy and interstitial dipole configurations and the 
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Figure 3. Flow stress for copper monocrystal: 1 - model data (solid line); 2 - experimental data (cross) 1211 . The numbers near curves are the 
deformation temperature given in Kelvin degrees. The strain rate is 10~ 4 i~' . 



concentration of interstitial atoms, vacancies, and divacancies were chosen to be zero. Fig.|3]shows strain hardening 
curves for copper monocrystal at various temperatures (solid lines) in comparison with experimental data El . A 
good agreement between experimental data and simulated curves is obtained. 

Fig. |4]2 shows strain hardening curves for copper monocrystal at various strain rate (solid lines) in in comparison 
with experimental data E2l . Fig. |4j? shows the square root of the dislocation density plotted versus the shear strain 
for copper monocrystal (solid line) in comparison with experimental data l23l . 
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Figure 4. Flow stress for copper monocrystal: model data (solid line) and experimental data (cross) 1221 . The numbers near curves are the strain 
rate given in s' 1 : 1-1.27- 10~ 2 , 2-4.18- 10" 3 , 3-4.34- 10" 4 , 4-4.09- 10" 5 , 5-4.46- 1(T 6 ; (b). The square root of the dislocation density 
plotted versus the shear strain for copper monocrystal: 1 - experimental data [23], 2 - model data. The initial dislocation density is l()' 2 m -2 and 
the strain rate is 10~ 4 .v~' . 

7. Conclusion 

In this work a mathematical model of plastic deformation in FCC metals was presented. Various features of plastic 
deformation by slip under various loading conditions can be described within this model. The model was implemented 
in-house software SPFCC (Slip Plasticity of Face-Centered Cubic). It provides graphic user interface and supports a 
study of plastic deformation of FCC metals. 

The "Slip Plasticity of Face-Centered Cubic" software creates a comfortable and friendly environment for compu- 
tational experiments while modelling plastic deformation processes in both FCC-metals and the dispersion-strengthened 
materials under various loading conditions. The proposed modelling technic automates many user actions, allowing to 
explore a variety of factors influencing plastic deformation of FCC metals and dispersion-strengthened materials. The 
results of the calculation are stored in the database of the software and are available for future reference. The software 
also includes a number of tools for collaboration within the scientific community. Tabular representation of the results 
of the calculations is designed to export data to third-party software for further post-processing. 

It should be noted that only the flow stress curves of copper with constant strain rate were considered in this work 
as an example. It was verified that the proposed mathematical model describes the plasticity of the FCC metals at 
various temperatures and strain rates accurately, and the results of computer simulation are in a good agreement with 
the experimental data. This computational investigation provides a fundamental understanding of the processes of 
intensity generation and annihilation of deformation defects of the centered cubic metals for constant strain rate. It's 
also demonstrates the influence of the typical parameters values of FCC metals on the plastic deformation features. 
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